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We consider spin dynamics for implementation in an atomistic framework and we address the 
feasibility of capturing processes in the femtosecond regime by inclusion of moment of inertia. In the 
spirit of an s-d-like interaction between the magnetization and electron spin, we derive a generalized 
equation of motion for the magnetization dynamics in the semi-classical limit, which is non-local in 
both space and time. Using this result we retain a generalized Landau-Lifshitz-Gilbert equation, 
also including the moment of inertia, and demonstrate how the exchange interaction, damping, and 
moment of inertia, all can be calculated from first principles. 
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In recent years there has been a huge increase in the in- 
terest in fast magnetization processes on a femto-second 
scale, which has been initiahzed by important develop- 
ments in experimental techniques [IHS], as well as po- 
tential technological applications [6]. From a theoretical 
side, the otherwise trustworthy spin dynamical (SD) sim- 
ulation method fails to treat this fast dynamics due to 
the short time and length scales involved. Attempts have 
been made to generalize the mesoscopic SD method to an 
atomistic SD, in which the dynamics of each individual 
atomic magnetic moment is treated [71 H]. While this 
approach should in principle be well suited to simulate 
the fast dynamics observed in experiments, it has not 
yet reached full predictive power as it has inherited phe- 
nomenological parameters, e.g. Gilbert damping, from 
the mesoscopic SD. The Gilbert damping parameter is 
well established in the latter regime but it is not totally 
clear how it should be transferred to the atomic regime. 
In addition, very recently it was pointed out that the mo- 
ment of inertia, which typically is neglected, plays an im- 
portant role for fast processes [9] . In this Letter we derive 
the foundations for an atomistic SD where all the rele- 
vant parameters, such as the exchange coupling, Gilbert 
damping, and moment of inertia, can be calculated from 
first principles electronic structure methods. 

Usually the spin dynamics is described by the phe- 
nomenological Landau-Lifshitz-Gilbert (LLG) equation 
[101 lllj which is composed of precessional and damping 
terms driving the dynamics to an equilibrium. By in- 
cluding the moment of inertia, we arrive at a generalized 
LLG equation 

M =M X (-7B + GM + iM) (1) 

where G and I are the Gilbert damping and the moment 
of inertia tensors, respectively. In this equation the effec- 
tive field B includes both the external and internal fields, 
of which the latter includes the exchange coupling and 
anisotropy effects. Here, we will for convenience include 
the anisotropy arising from the classical dipole-dipole in- 
teraction responsible for the shape anisotropy as a part 



of the external field. The damping term in the LLG 
equation usually consists of a single damping parame- 
ter, which essentially means that the time scales of the 
magnetization variables and the environmental variables 
are well separated. This separation naturally brings a 
limitation to the LLG equation concerning its time scale 
which is restricting it to the mesoscopic regime. 

The addition of a moment of inertia term to the LLG 
equation can be justifies as follows. A general process 
of a moment M under the influence of a field F is al- 
ways endowed with inertial effects at higher frequencies 
|12) . The field F and moment M can, for example, be 
stress and strain for mechanical relaxation, electric field 
and electric dipole moment in the case of dielectric re- 
laxation, or magnetic field and magnetic moment in the 
case of magnetic relaxation. In this Letter we focus on 
the latter case — the origin of the moment of inertia in 
SD. The moment of inertia leads to nutations of the mag- 
netic moments, see Fig. [l] Its wobbling variation of the 
azimuthal angle has a crucial role in fast SD, such as fast 
magnetization reversal processes. 

In the case of dielectric relaxation the inertial effects 
are quite thoroughly mentioned in the literature [121 [H] , 
especially in the case of ferroelectric relaxors. Coffey et 
al. [HI have proposed inertia corrected Debye's theory of 
dielectric relaxation and showed that by including inertial 
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FIG. 1: The three contributions in Eq. ([T|, the bare preces- 
sion arising from the effective magnetic field, and the super- 
imposed effects from the Gilbert damping and the moment of 
inertia, respectively. 
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effects, the unphysical high frequency divergence of the 
absorption co-efficient is removed. 

Very recently Ciornei et al [5] have extended the LLG 
equation to include the inertial effects through a mag- 
netic retardation term in addition to precessional and 
damping terms. They considered a collection of uni- 
formly magnetized particles and treated the total angular 
momentum L as faster variable. They obtained Eq. ([T]) 
from a Fokker-Plank equation where the number den- 
sity of magnetized particles were calculated by integrat- 
ing a non-equilibrium distribution function over faster 
variables such that faster degrees of freedom appear as 
parameter in the calculation. 

The authors showed that at very short time scales the 
inertial effects become important as the precessional mo- 
tion of magnetic moment gets superimposed with nuta- 
tion loops due to inertial effects. It is pointed out that 
the existence of inertia driven magnetization dynamics 
open up a pathway for ultrafast magnetic switching [15 
beyond the limitation 16J of the precessional switching. 

In practice, to perform atomistic spin dynamics simu- 
lations the knowledge of G and I is necessary. There are 
recent proposals [T71 [TH] of how to calculate the Gilbert 
damping factor from first principles in terms of Kubo- 
Greenwood like formulas. Here, we show that similar 
techniques may by employed to calculate the moment of 
inertia tensor I. Finally, we present a microscopical jus- 
tification of Eq. ([T]), considering a collective magnetiza- 
tion density interacting locally with electrons constitut- 
ing spin moments. Such a description would in principle 
be consistent with the study of magnetization dynam- 
ics where the exchange parameters are extracted from 
first-principles electronic structure calculations, e.g den- 
sity functional theory (DFT) methods. We find that in 
an atomistic limit Eq. ([T]) actually has to be general- 
ized slightly as both the damping and inertia tensors are 
naturally non-local in the same way as the exchange cou- 
pling included in the effective magnetic field B. From 
our study it is clear that both the damping and the mo- 
ment of inertia effects naturally arise from the retarded 
exchange interaction. 

We begin by considering the magnetic energy _E = M • 
B. Using that its time derivative isi; = M- B-|-M-B 
along with Eq. ([T]), we write 

E =M • B + • [ciA + IM^ . (2) 

Relating the rate of change of the total energy to the 
Hamiltonian "H, through E = [dT-L/dt), and expanding 
the % linearly around its static magnetization Mq, with 
M.{t) = Mo we can write nK.UQ + fi{t) ■ V^-Hq, 

where Hq = 'H(Mo). Then the rate of change of the total 
energy equals E = fi - (V^H) to the first order. Following 
Ref. [12] and assuming sufficiently slow dynamics such 
that ^lit') = ^(t) — Tfi{t) + T^il{t)/2, T = t — t', wc can 



write the rate of change of the magnetic energy as 

E = hm^/ij[xij(w)^j + id^Xt]{^)i^] - dlxiji^Yl^] /'A- 

(3) 

Here, x^J{^) = /(-z)0(r)([a,Ho(t), a,Ho(t')])e'"^di', 
r = t — t' , is the (generalized) exchange interaction 
tensor out of which the damping and moments of in- 
ertia can be extracted. Summation over repeated in- 
dices = x,y,z) is assumed in Eq. ([3|. Equat- 
ing Eqs. ([2]) and ([s]) results in an internal contribu- 
tion to the effective field about which the magnetiza- 
tion precesses Bi^t = ^HuIi^^q xiuj), the damping term 
G = 7 limi^_j.o «9a;X('^) a-s well as the moment of inertia 
i = -7lim^^oa2;^(a;)/2. 

For a simple order of magnitude estimate of the damp- 
ing and inertial coefficients, G and I, respectively, we 
may assume for a state close to a ferromagnetic state 
that the spin resolved density of electron states ^^(e) 
corresponding to the static magnetization configuration 
Hq is slowly varying with energy. At low temperatures 
we, then, find 

G ~277rsp [{d,Ho)p{djHo)p]^^^^ , (4) 

in agreement with previous results [l9j. Here, sp denotes 
the trace over spin 1/2 space. By the same token, the 
moment of inertia is estimated as 

i ^ - (j/D) sp mHo)p{d,Ho)pl^^^ , (5) 

where 2D is the band width of density of electron states 
of the host material. Typically, for metallic systems the 
band width 2D ~ 1 — 10 eV, which sets the time-scale 
of the inertial contribution to the femto second (10~^^ 
s) regime. It, therefore, defines magnetization dynamics 
on a time-scale that is one or more orders of magnitude 
shorter compared to e.g. the precessional dynamics of the 
magnetic moment. 

Next, we consider the physics leading to the LLG equa- 
tion given in Eq. ([!]) . As there is hardly any microscopical 
derivation of the LLG equation in the literature, we in- 
clude here, for completeness the arguments that leads to 
the equation for the spin-dynamics from a quantum field 
theory perspective. 

In the atomic limit the spin degrees of freedom are 
deeply intertwined with the electronic degrees of free- 
dom, and hence the main environmental coupling is the 
one to the electrons. In this study we are mainly con- 
cerned with a mean field description of the electron 
structure, as in the spirit of the DFT. Then a natural 
and quite general description of the magnetic interac- 
tion due to electron-electron interactions on the atomic 
site around r within the material is captured by the s-d- 
like model "Hint = — J ^(r, r')M(r, t) ■ s(r', t)drdr' , where 
J(r,r') represents the interaction between the magneti- 
zation density M and the electron spin s. From a DFT 



3 



perspective the interaction parameter J(r,r') is related 
to the effective spin dependent exchange-correlation func- 
tional Bxc[M{r')]{r). For generality we assume a fully 
relativistic treatment of the electrons, i.e. including the 
spin-orbit coupling. In this interaction the dichotomy of 
the electrons is displayed, they both form the magnetic 
moments and provide the interaction among them. 

Owing to the general non-equilibrium conditions in the 
system, we define the action variable 



moment is kept constant. We, thus, find 
M(r,i) =M(r,i) x ( --fB,^t{r,t) 



S — (j) Hint dt + Sz + 5wZWN 



(6) 



on the Keldysh contour [20H22j . Here, the ac- 
tion Sz = — 7/p/Bext(r,0 • 'M{r,t)dtdr represents 
the Zeeman coupling to the external field Boxt(r,i), 
whereas the Wess-Zumino-Witten-Novikov (WZWN) 
term Swzwn = / /g^ M(r, r) • [drM{r,t;T) x 
9tM(r, t; r)](iTc?t|M(r)|~^(ir describes the Berry phase 
accumulated by the magnetization. 

In order to acquire an effective model for the magne- 
tization density M(r, t), we make a second order [27 ex- 
pansion of the partition function Z[M(r, t)] = tr Tce^'^ , 
and take the partial trace over the electronic degrees of 
freedom in the action variable. The effective action SS]\{ 
for the magnetization dynamics arising from the mag- 
netic interactions described in terms of "Hint, can, thus, 
be written 

6Sm=- <j> J M{r, t) ■V{r,r'- t,t') -Mir', t')drdr'dtdt', 

(7) 

where V{r,r';t,t') = / J(r, ri)(-i)(Ts(ri, t)s(r2, t')> x 
J(r2, r')dridr2 is a dyadic which describes the electron 
mediated exchange interaction. 

Conversion of the Keldysh contour integrations into 
real time integrals on the interval (— oo, oo) results in 

J M(f''"*)(r,i) • [M(r,t) x M{r,t)]dt\M{r)\~^dr 

+ /M<-..,..,^I.',..';UrM(r',n**'<i«*' 

-J I Bext (r , t) • M(f-^') (r , t)dtdr, (8) 

with M^f'^"') (r, t) = Mu{r,t) - M/(r,i) and M{r,t) = 
[M„(r,i) -I- M;(r,t)]/2 which define fast and slow vari- 
ables, respectively. Here, M^(;) is the magnetization den- 
sity defined on the upper (lower) branch of the Keldysh 
contour. Notice that upon conversion into the real time 
domain, the contour ordered propagator V is replaced by 
its retarded counterpart V^. 

We obtain the equation of motion for the (slow) mag- 
netization variable ]VI(r,t) in the classical limit by mini- 
mizing the action with respect to M(^'^''*)(r,i), cross mul- 
tiplying by 'M.{r,t) under the assumption that the total 



J V''{r,r';t,t') ■M{r',t')dt'dr' 



(9) 



Eq. ([9]) provides a generalized description of the semi- 
classical magnetization dynamics compared to the LLG 
Eq. ([1]) in the sense that it is non-local in both time and 
space. The dynamics of the magnetization at some point 
r depends not only on the magnetization locally at r, 
but also in a non-trivial way on the surrounding magne- 
tization. The coupling of the magnetization at different 
positions in space is mediated via the electrons in the 
host material. Moreover, the magnetization dynamics is, 
in general, a truly non-adiabatic process in which the 
information about the past is crucial. 

However, in order to make connection to the magne- 
tization dynamics as described by e.g. the LLG equa- 
tion as well as Eq. ([T]) above, we make the following 
consideration. Assuming that the magnetization dy- 
namics is slow compared to the electronic processes in- 
volved in the time-non-local field T>{r,r';t,t'), we ex- 
pand the magnetization in time according to M(r',t') w 
M(r',t) - TM(r',i) + T'^M{r' ,t)/2. Then for the inte- 
grand in Eq. we get 

P'-(r,r';i,<')-M(r',i') = 

V^{r, r'; t, t') ■ [M(r', t) - rM(r', t) + y M(r', t)]. (10) 

Here, we observe that as the exchange coupling for the 
magnetization is non-local and mediated through T), this 
is also true for the damping (second term) and the inertia 
(third term). 

In order to obtain an equation of the exact same 
form as LLG in Eq. ([!]) we further have to assume 
that the magnetization is close to a uniform ferromag- 
netic state, then we can justify the approximations 
M(r',t) w M(r,t) and M(r',i) w M(r,t). When 
Bint = ~ SV{Y,Y'-t,t') ■ M{r', t)dr'dt'/"f is included in 
the total effective magnetic field B, the tensors of Eq. ([T]) 
G and I can be identified with — / T'D{r,r';t,t')dr'dt' 
and / T'^'D{r,r';t,t')dr'dt' /2, respectively. From a first 
principles model of the host materials we have, thus, de- 
rived the equation for the magnetization dynamics dis- 
cussed in Ref. where it was considered from purely 
classical grounds. However it is clear that for a treatment 
of atomistic SD that allows for all kinds of magnetic or- 
ders, not only ferromagnetic, Eq. ([ij is not sufficient and 
the more general LLG equation of Eq. (|9| together with 
Eq. ( 10 ) has to be used. 

We finally describe how the parameters of Eq. ([T]) 
can be calculated from a first principles point of view. 
Within the conditions defined by the DFT system, the 
interaction tensor T>^ is time local which allows us to 
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write lime_^o«42?''(r,r';£) = J {r,r';t,t')dt' and 
lime^oc'e^''(r,r';£) = - / T^V{v,v']t,t')dt' , where 



I?'-(r,r';e) =4sp / Jr^Jp, 



e — cj - 



^(5 



X ^ImGPp(c.)^ImG^^, ^ ^dpdp'. (11) 



27r 27r 



Here, Jrr' = J(r, r') whereas G^^,{uj) = G^{r,r';uj) is 
the retarded GF, represented as a 2 x 2-matrix in spin- 
spaces. We notice that the above result presents a general 
expression for frequency dependent exchange interaction. 
Using Kramers-Kronig's relations in the limit e — )• 0, it 



is easy to see that Eq. ( 11 ) leads to 
2?'-(r,r';0)= ^ 



Splm J JrpJp'r'f{<^) 

X cTGl,Auj)cTGlA^)dujdpdp', (12) 



in agreement with e.g. Ref. ^24j. We can make connection 
with previous results, e.g. Refs. [13 UHl and observe that 



Eq. (11) contains the isotropic Heisenberg, anisotropic 



Ising, and Dzyaloshinsky-Moriya exchange interactions 
between the magnetization densities at different points 
in space [H], as well as the onsitc contribution to the 
magnetic anisotropy. 



Using the result in Eq. ( 11 ), we find that the damping 



tensor is naturally non-local and can be reduced to 



G(r,r') = 



1 



JrpJp'r' f' {^) 



cTlmG'p,Joj)<TlinGl ,{Lj)dLjdpdp' , (13) 



which besides the non-locality is in good accordance with 
the results in Refs. [T7| US], and is closely connected to 
the so-called torque-torque correlation model [57]. With 
inclusion of the the spin-orbit coupling in C", it has been 



demonstrated that Eq. ( 13 1 leads to a local Gilbert damp- 



ing of the correct order of magnitude for the case of fer- 
romagnetic permalloys [17i . 

Another application of Kramers-Kronig's relations 
leads, after some algebra, to the moment of inertia tensor 



I(r,r')=sp J JrpJp'r'fHrT[lraG-p,piu:)adiRcG-pp,{oj] 

+ ImG^p, (c.)^a2 j^eGPp(w)] ^dpdp\ (14) 

where we notice that the moment of inertia is not sim- 
ply a Fermi surface effect but depends on the electronic 
structure as a whole of the host material. Although the 
structure of this expression is in line with the exchange 



coupling in Eq. (12 1 and the damping of Eq. (13 1, it is 



a little more cumbersome to compute due the presence 
of the derivatives of the Green's functions. Note that it 
is not possible to get completely rid of the derivatives 
through partial integration. These derivatives also make 



the moment of inertia very sensitive to details of the elec- 
tronic structure, which has a few implications. Firstly the 
moment of inertia can take large values for narrow band 
magnetic materials, such as strongly correlated electron 
systems, where these derivatives are substantial. For 
such systems the action of moment of inertia can be im- 
portant for longer time scales too, as indicated by Eq. ([5|. 
Secondly, the moment of inertia may be strongly depen- 
dent on the reference magnetic ordering for which it is 
calculated. It is well known that already the exchange 
tensor parameters may depend on the magnetic order. 
It is the task of future studies to determine how trans- 
ferable the moment of inertia tensor as well as damping 
tensor are in-between different magnetic ordering. 

In conclusion, we have derived a method for atomistic 
spin dynamics which would be applicable for ultrafast 
(femtosecond) processes. Using a general s-d-like interac- 
tion between the magnetization density and electron spin, 
we show that magnetization couples to the surrounding 
in a non-adiabatic fashion, something which will allow for 
studies of general magnetic orders on an atomistic level, 
not only ferromagnetic. By showing that our method 
capture previous formulas for the exchange interaction 
and damping tensor parameter, we also derive a formula 
for calculating the moment of inertia from first principles. 
In addition our results point out that all parameters are 
non-local as they enter naturally as bilinear sums in the 
same fashion as the well established exchange coupling. 
Our results are straight-forward to implement in existing 
atomistic SD codes, so we look on with anticipation to 
the first applications of the presented theory which would 
be fully parameter-free and hence can take a large step 
towards simulations with predictive capacity. 
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